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Abstract 

A new universal theory for flow instability and turbulent transition is proposed in 
this study. Flow instability and turbulence transition have been challenging subjects for 
fluid dynamics for a century. The critical condition of turbulent transition from theory 
and experiments differs largely from each other for Poiseuille flows. This enigma has 
not been clarified so far owing to the difficulty of the problem. In this paper, a new 
mechanism of flow instability and turbulence transition is presented for parallel shear 
flows and the energy gradient theory of hydrodynamic instability is proposed. It is stated 
that the total energy gradient in the transverse direction and that in the streamwise 
direction of the main flow dominate the disturbance amplification or decay. Thus, they 
determine the critical condition of instability initiation and flow transition under given 
initial disturbance. A new dimensionless parameter K for characterizing flow instability 
is proposed for wall bounded shear flows, which is expressed as the ratio of the energy 
gradients in the two directions. It is thought that flow instability should flrst occur at 
the position of Krnax which may be the most dangerous position. This speculation is 
confirmed by Nishioka et al's experimental data. Comparison with experimental data 
for plane Poiseuille flow and pipe Poiseuille flow indicates that the proposed idea is really 
valid. It is found that the turbulence transition takes place at a critical value of -fCmax 
of about 385 for both plane Poiseuille flow and pipe Poiseuille flow, below which no 
turbulence will occur regardless the disturbance. More studies show that the theory is 
also valid for plane Couette flows and Taylor-Couette flows between concentric rotating 
cylinders. It is concluded that the energy gradient theory is a universal theory for the 
flow instability and turbulent transition which is valid for both pressure and shear driven 
flows in both parallel flow and rotating flow configurations. 
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1 Introduction 



Understanding the mechanism of turbulence has been a great challenge for over a century. 
Now, it is still very far from approaching a comprehensive theory and the final resolution of 
the turbulent problem [1]. Reynolds (1883) [1] [2] did the first famous experiments on pipe 
flow demonstrating the transition from laminar to turbulent flows. Since then, various stabil- 
ity theories emerged during the past 120 years for this phenomenon, but few are satisfactory 
in the explanation of the various flow instabilities and the related complex flow phenomena. 

The pipe Poiseuille flow (Hagen-Poiseuille) is linearly stable for all the Reynolds number 
Re by eigenvalue analysis [3] [4] [5] [6] [7]. However, experiments showed that the flow would 
become turbulence if Re (= pUD/ fi) exceeds a value of 2000. Experiments also showed that 
disturbances in a laminar flow could be carefully avoided or considerably reduced, the onset 
of turbulence was delayed to Reynolds numbers up to Re = O(IO^) [6] [7]. For Re above 
2000, the characteristic of turbulence transition depends on the disturbance amplitude and 
frequency, below which transition from laminar to turbulent state does not occur regardless 
of the initial disturbance amplitude [8] [9]. Thus, it is clear that the transition from laminar 
to turbulence for Re > 2000 is dominated by the behaviours of mean flow and disturbance. 
Only the combined effect of the two factors reaches the critical condition, could the transition 
occur for Re > 2000. These are summarized in Tabic 1. 

Linear stability analysis of plane parallel flows gives critical Reynolds number Re (= 
pVoh/ jj) of 5772 for plane Poiseuille flow, while experiments show that transition to turbu- 
lence occurs at Reynolds number of order 1000 [10] [3] [6] [7], even though the laminar flow 
could also kept to Re = O(IO^) [11]. One resolution of these paradoxes is that the domain of 
attraction for the laminar state shrinks for large Re (as Re'*' say, with 7 < 0), so that small 
but finite perturbations lead to transition [6] [12]. Grossmann [7] commented that this dis- 
crepancy demonstrates that nature of the onset-of-turbulence mechanism in this flow must 
be different from an eigenvalue instability. Orszag and Patera [13] remarked that the mecha- 
nism of transition is not properly represented by parallel-flow linear stability analysis. They 
proposed a linear three-dimensional mechanism to predict the transitional Reynolds number. 
Some nonlinear stability theories have been proposed, for example, in [14]. However, these 
theories do not seem to offer a good agreement with the experimental data. 

Energy method was also used in the study of flow instabilities [15] [16] [17] [4] [18] [5]. 
In energy method, one observes the rate of increasing of disturbance energy to study the 
instability of the flow system. The critical condition is determined by the maximum Reynolds 
number at which the disturbance energy in the system monotonically decreases. In the flow 
system, it is considered that turbulence shear stress interacts with the velocity gradient 
and the disturbance gets energy from mean flow in such a way. Thus, the disturbance is 
amplified and the instability occurs with the energy increasing of disturbance. Therefore, it is 
recognized that it is the basic state vorticity leading to instability. The energy method could 
not get agreement with the experiments either [16] [4] [5]. In recent years, various transition 
scenarios have been proposed [6] [7] [19] [20] [21] [22] [12] [23] for the subcritical transition. 
Although we can get a better understanding of the transition process from these scenarios, 
the mechanism is still not fully understood and the agreement with the experimental data 
is still not satisfied. 

Generally, the transition from laminar flow to turbulent flow is not generated suddenly 
in the entire flow field but it first starts from somewhere in the flow field and then spreads 
out gradually from this position. As is well known in solid mechanics, the damage of a 
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Base flow 


Disturbance 


Resulting flow state 


Re < 2000 


No matter how large. 


Disturbance decay; flow keeps laminar. 


Re ~ 10^ 


Kept low. 


Flow keeps laminar. 


Re > 2000 


Enough large, depend on Re. 


Transition occurs. 



Table 1: Characteristic of transition of Hagen-Poiseuille flow (pipe Poiseuille flow). 



metal component generally starts from some area such as manufacturing fault, crack, stress 
concentration, or fatigue position, etc. In fluid mechanics, we consider that the breaking 
down of a steady laminar flow should also start from a most dangerous position first. The 
consequent questions are: (a) Where is this most dangerous position for Poiseuille flow? (b) 
What is the mechanism and the dominant factor for this phenomenon? (c) What parameter 
should he used to characterize this position? These questions are our concern. Finding the 
solution of these problems is important to the understanding of the phenomenon and the 
estimation of flow transition. Because the turbulence transition is generally resulted in by 
flow instability [15], we think that the critical condition of transition should be determined by 
the position where the flow instability first takes place. If the mechanism of flow instability 
is sought out and the most dangerous position is found, the critical condition of transition 
could be determined. 

In this study, wc explore the critical condition of main flow for instability and turbulence 
transition, and not deal with the detailed process of disturbance amplification. The energy 
gradient theory is proposed to explain the mechanism of flow instability and turbulence 
transition for parallel flows. A new dimensionless parameter for characterizing the critical 
condition of flow instability is proposed. Comparison with experimental data for plane 
Poiseuille flow and pipe Poiseuille flow at subcritical transition indicates that the proposed 
idea is really valid. 

2 Proposed Mechanism of Flow Instability 

The plane Poiseuille flow in a channel is shown in Fig.l. For the given flow geometry and 
fluid property, with the increasing mean velocity [/, the flow may transit to turbulence if 
the Re exceeds a critical value (under certain disturbance). The velocity profiles for laminar 
and turbulent flows are shown respectively in Fig. 2. It can be imagined that there is a 
"driving factor" which pulls the laminar velocity profile outward toward the walls when 
the transition takes place. What should be such a driving factor? From a large amount 
of observations, it is thought that the increase of the gradient of fiuid kinematic energy in 
transverse direction, ^(|^^) = form a driving force" to cause the increase of 

flow disturbance for given flow condition, while the gradient of the viscous friction force may 
resist or absorb the disturbance. Here, V is the magnitude of local velocity. The stability of 
the flow depends upon the effects of these two roles. With the increasing of mean velocity 
U for parallel flows, the energy gradient in the transverse direction increases. If this energy 
gradient is large enough it will lead to a disturbance ampliflcation of the flow. The viscosity 
friction caused by shear stress would stabilize the flow by absorbing the velocity fluctuation. 
When the energy gradient in transverse direction reaches beyond a critical value, the laminar 
flow could not balance this disturbance and flow instability might be excited. Finally, the 
turbulence flow would be triggered when the transverse energy gradient continuously keeps 
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Figure 1: Velocity distribution variation with the increased Reynolds number for given fluid 
and geometry in plane Poiseuille flows. Re = pUL/ ^x, L = 2h, where h is the half- width of 
the channel. 




large enough with the flow forward. The energy gradient in the transverse direction makes 
the exchange of energy between the fluid layers and sustains the turbulence. Therefore, it is 
proposed that the necessary condition for the turbulence transition is that there is an energy 
gradient in the transverse direction of the main flow. 

Now, we prove this necessary condition is correct at least for parallel flows. If the gravita- 
tional energy is neglected, the total energy gradient in transverse direction is ^ (p + ^pV^) . 

For parallel flows, || = and V = u. If this energy gradient is zero, ^(5^^) = = 0) 
there must be ^ = due toV^O. Thus, the rate of increase of disturbance energy will be 
negative due to viscous dissipation because the disturbance could not obtain energy from the 
base flow owing to zero velocity gradient [16] [4] [5]. Therefore, the disturbance must decay 
at this case. In such way, it is proved that the energy gradient in the transverse direction is 
a necessary condition for the flow transition. 

In addition, when there is a pressure gradient in the normal direction to the flow direction, 
this pressure gradient could also result in a flow instability even the Reynolds number is low. 
Both centrifugal and Coriolis instabilities are those caused by pressure gradients. Elastic 
instability is also that produced by the transversal pressure gradient [25] [26] . The mechanism 
of instability should take into account of the effect of the variation of cross-streamline pressure 
for these cases, which may lead to flow instability or accelerates the instability initiation. In 
some cases of incompressible flows such as stratified flows, the gravitational energy should 
be taken into account. 

For given flow geometry and fluid, it is proposed that the flow stability condition can be 
expressed as, 

^(^pgyy+P+\pV^^ <C, (1) 

where gy is the component of gravity acceleration in y direction, and C is a constant which 
is related to fluid property and geometry. The x axis is along the flow direction and the y 
axis is along the transverse direction. In this study, we first show that the proposed idea is 
really correct for Poiseuille flows. 

3 Formulation and Theory Description 

The conservation of momentum for an incompressible Newtonian fluid is (neglecting gravity 
force): 

^'^~dt + V • VV) = -Vp + /xV^ V. (2) 
V • VV = (V • V) - V X V X V, (3) 



Using the identity. 



equation (2) can be rearranged as, 

p^ + Vip+^pV^)=pV^\ + p(VxVx\), (4) 

where p is the fluid density, t the time, V the velocity vector, p the hydrodynamic pressure, 
p the dynamic viscosity of the fluid. If the viscous force is zero, the above equation becomes 
the Lamb form of momentum equation. This equation can be found in most text books. For 
incompressible flow, the total pressure represents the total energy in Eq. (4). Actually, the 



5 



energy equation has long been used for stability analysis as previously mentioned [18] [4] [5] 
[17] [16]. 

In previous sections, it is proposed that the instability of viscous flows depends on the 
relative magnitude of the energy gradient in transverse direction and the viscous friction 
term. A larger energy gradient in transverse direction tries to lead to amplification of a 
disturbance, and a large shear stress gradient in streamwise direction tends to absorb this 
disturbance and to keep the original laminar flow state. The transition of turbulence depends 
on the relative magnitude of the two roles of energy gradient amplification and viscous friction 
damping under given disturbance. We propose the parameter for characterizing the relative 
role of these effects below. 

Let ds represent the differential length along a streamline in a Cartesian coordinate 
system, 

ds = dx + dy. (5) 



With dot multiplying Eq. (4) by ds, we obtain, 

p— ■ ds + V(p + -pV^) ■ ds = pV^Y • ds + p(V X V X V) • ds. (6) 



Since (V x V x V) • ds =0 along the streamline, for steady flows, we obtain the energy 
gradient along the streamline, 

1 rja 

dip + i^pv'')/ds = pV'V ■ —={pV'y)^. (7) 

This equation shows that the total energy gradient along the streamwise direction equals to 
the viscous term (//V^V)^. For pressure driving flows, (/xV^V)^ represents the energy loss 
due to friction. It is obvious that the total energy decreases along the streamwise direction 
due to viscous friction loss. The energy gradient along the transverse direction is, 

a(p + i,vVa„ = | + pKf. (8) 

It can be seen that the energy gradient at transverse direction depends on the velocity 
gradient and the velocity magnitude as well as the transversal pressure gradient. 

The relative magnitude of the energy gradients in the two directions can be expressed 
by a new dimensionless parameter, K, the ratio of the energy gradient in the transverse 
direction to that in the streamwise direction, 

_ dE/dn _ d{p + yv^)/dn _ dp/dn + pV{dV/dn) 

~ dE/ds ~ dip + \pV^)/ds ~ (/xV^V), ' ^' 

where E = p + ^pV^ expresses the total energy, n denotes the direction normal to the 
streamwise direction, and s denotes the streamwise direction. 

It is noticed that the parameter K is a field variable and it represents the direction of 
the vector of local total energy gradient. It can be seen that when K is small, the role of the 
viscous term in Eq.(9) is large and the flow tends to be stable. When K is large, the role of 
numerator in Eq.(9) is large and the flow tends to be unstable. For given flow fleld, there 
is a maximum of K in the domain which represents the most possible unstable location. 



6 



Therefore, in the area of high value of the flow tends to be more unstable than that in 
the area of low value of K. The first instability should be initiated by the maximum of K, 
Kmaxi in the flow field for given disturbance. In other words, the position of maximum of 
K is the most dangerous position. The magnitude of Kmax in the flow field symbolizes the 
extent of the flow approaching the instability. Especially, if Kmax = oo, the flow is certainly 
potential to be unstable under some disturbance. For given flow disturbance, there is a 
critical value of Kmax over which the flow becomes unstable. In particular, corresponding to 
the subcritical transition in wall bounded parallel flows, the Kmax reaches its critical value, 
below which no transition occurs regardless of the disturbance amplitude. For unidirectional 
parallel flows, this critical value of Kmax should be a constant regardless of the fluid property 
and the magnitude of the geometrical parameter. Now, owing to the complexity of the flow, 
it is difficult to predict this critical value by theory. Nevertheless, it can be determined by 
the experimental data for given flows. 

For parallel flows, the coordinates s — n becomes the x — y coordinates and the pressure 
gradient dp/dn = dp/dy = 0. Thus, using the global approximation, d{^pV^)/dy ^ pU^/L, 
(//V^V)^ - d'^u/dy^ - U/L^, we obtain K ^ pUL/fi = Re, where L is a characteristic 
length in the transverse direction. Thus, the parameter K is equivalent to the Reynolds 
number in the global sense for parallel flows. 

The negative energy gradient in the streamwise direction plays a part of resisting on or 
absorbing the disturbance. Therefore, it has a stable role to the flow. The energy gradient 
in the transverse direction has a role of amplifying disturbance. Therefore, it has a unstable 
role to the flow. Here, it is not to say that a high energy gradient in transverse direction 
necessarily leads to instability, but it has a potential for instability. Whether instability 
occurs also depends on the magnitude of the disturbance. From this discussion, it is easily 
understood that the disturbance amplitude required for instability is small for high energy 
gradient in transverse direction (high Re) [9]. 

The value of "arctani^" expresses the angle between the direction of the total energy 
gradient and the streamwise direction. We write that, 

a = arctanif. (10) 

In this paper, the angle a is named as "energy angle," as shown in Fig.3. The value of the 
energy angle (its absolute value) can also be used to express the extent of the flow near 
the instability occurrence. Fig.3 and Fig. 4 show the schematic of the energy angle for some 
flows. There is a critical value of energy angle, ac, corresponding to the critical value of 
Kmax- When a > ac, the flow becomes unstable. For Poiseuille flows, 0° < a < 90° and the 
stability depends on the magnitudes of the energy angle and the disturbance. For parallel 
flow with a velocity inflection, a = 90° (Kmax = oo) at the inflection point and the flow is 
therefore unstable. This is for the first time to theoretically prove that viscous flow with a 
velocity inflection is unstable. Inversely, as discussed before, viscous flow without velocity 
inflection may not be necessarily stable, depending on the flow conditions (as shown in table 
1). These results are correct at least for pressure driving flows. 

Fig.4 is a best explanation of inflectional instability for viscous flows. According to 
Eq.(9), the physics of the criterion presented in this paper is easily understood. The energy 
gradient in the transverse direction tries to amplify a small disturbance, while the energy 
loss due to friction in the streamwise direction plays a damping part to the disturbance. The 
parameter K represents the relative magnitude of disturbance amplification due to energy 
gradient and disturbance damping of viscous loss. When there is no inflection point in the 
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0° < a < 9Q 




Figure 4: Schematic of the direction of the total energy gradient and energy angle for flow 
with an inflection point at which the energy angle equals 90 degree. 
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velocity distribution, the amplification or decay of disturbance depends on the two roles 
above, i.e., the value of K. When there is an inflection point in the velocity distribution, the 
viscous term vanishes, and while the transversal energy gradient still exists at the position 
of inflection point. Thus, even a small disturbance must be amplified by the energy gradient 
at this location. Therefore, the flow will be unstable. Rayleigh (1880) [24] only proved that 
inviscid flow with an inflection point is unstable, while wc demonstrate here that viscoTis flow 
with an inflection point is unstable. The Raylcigh's criterion was derived from mathematics 
using the inviscid theory, but its physical meaning is still not clear. 

It is well known that viscosity plays dual role on the flow instability and disturbance 
amplification [4] [18] [5] [15]. It can be seen from equation (9) that the higher the viscosity, 
the larger the viscous friction loss. Thus, the flow is more stable. If the values of the vorticity 
and the streamwise velocity are high, the transversal energy gradient will be high. Thus, 
the flow is more unstable. Prom these discussion, it is known that viscosity mainly plays a 
stable role to the initiation of flow instability at subcritical transition by affecting the base 
flow through the viscous friction of streamwise velocity. This is consistent with criterion of 
the Reynolds number. 

Reynolds number represents the ratio of convective inertia force to the viscous force in 
the Navier-Stokes equations as a dimensionless parameter. However, the magnitude of the 
Reynolds number is only a global indication of the flow states and a rough expression for the 
transition condition. At same Re number, the behaviour of flow instability may be different 
due to the different combination of the magnitude of viscosity with other parameters when 
Re is larger than an indifference Reynolds number, as shown by the chart of the solution 
of the Orr-Sommcrfeld equation [18] [27] [28]. The role of viscosity is complicated with 
the variations of flow parameters and flow conditions. The generation of turbulence is not 
simply caused by increasing the Re, but actually by increasing the value of K. When Re 
increases, it inevitably leads to the increase of K in the flow field. The magnitude of local 
disturbance for a fixed point depends not only on the apparent Re number, but also on the 
local fiow conditions. For steady Poiseuille flow, the convective inertia term is zero and the 
flow becomes turbulent at high Re. We can see that the occurrence of turbulence is not due 
to the convective inertia term in this case. In Poiseuille flow, the local Reynolds number is 
high in the core area along the centerline, but the degree of turbulence is low. In the area 
near the wall, the local Reynolds number is low, but the degree of turbulence is high. In 
uniform flows, high Reynolds number may not necessarily leads to turbulence. In summary, 
the Reynolds number is only a global parameter, and the K is a local parameter which 
represents the local behaviour of the flow and best reflects the role of energy gradient. 



4 Analysis on Poiseuille Flows 
4.1 Plane Poiseuille Flow 

The instability in Poiseuille flows with increasing mean velocity U can be demonstrated as 
below for a given fliiid and flow geometry by using the energy gradient concept. 
For 2D Poiseuille flow, the momentum equation is written as, 

dp d'^u , . 

showing that viscous force term is proportional to the streamwise pressure gradient. The 
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Figure 5: Kinematic energy gradient and viscous force term versus the axial pressure gradient 
for plane Poiseuille flow for any position in the flow field. \dpldx\ oc [/ oc i?e for given fluid 
and geometry. 



integration to above equation gives, 



u = Uq\1 — 



where uq = —^h'^^i is the center line velocity, and h is the channel half width. 

Although the energy gradient is not explicitly shown in above equation, it can be ex- 
pressed as below for any position in the flow field (noticing v = 0), 



Thus, for any position in the flow field, 



1 



r 



dp 
dx 



dy \2 



oc 



dp 
dx 



(13) 



(14) 



The behaviours of the equations (11-14) are shown in Fig.5: the viscous force term is 

linear, while the kinematic energy gradient increases quadratically with the pressure gradient. 
Therefore, at low value of mean velocity U , the viscous friction term could balance the 
disturbance amplification caused by energy gradient. At large value of U , the viscous friction 
term may not constrain the disturbance amplification caused by energy gradient and the flow 
may transit to turbulence. 

For plane Poiseuille flows, the ratio of the energy gradient to the viscous force term, K, 
is [dp/dy = 0), 



K = 



d_ 
dy 



1 



d'^u 
dy^ 



'IpyuQ 



-Uo 



-fi- 



2no 
/i2 
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Figure 6: Velocity, energy, and K along the transversal direction y/h for plane Poiseuille 
flow, which are normalized by the their maximum. 



velocity at centerline and U is the averaged velocity. It can be seen that K is proportional 
to Re for a fixed point in the flow field. 

The distribution of u, E, and K along the transversal direction for plane Poiseuille flow 
is shown in Fig. 6. It is clear that there are maximum of K at y/h = ±0.5774, as shown in 
Fig. 6. This maximum can also be obtained by differentiating the equation (15) with y/h 
and letting the derivatives equal to zero. Since we concern the maginitude of the K and the 
velocity profile is symmetrical to the centerline, we refer the maximum of K as its positive 
value thereafter. We think that the flow breakdown of the Poiseuille flow should not suddenly 
occur in the entire flow field, but it first takes place at the location of i^max in the domain, 
then it spreads out according to the distribution of K value. The formation of turbulence 
spot in shear flows may be related to this procedure. 

4.2 Pipe Poiseuille Flow 

Similar analysis to the plane Poiseuille flow can be carried out for the circular Poiseuille flow. 
For circular pipe Poiseuille flow, the momentum equation is written as. 




(15) 



Here, Re = pUL/jj, and U = |iio has been used for plane Poiseuille flow, uq is the maximum 




(16) 
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Figure 7: Velocity, energy, and K along the transversal direction r/R for pipe Poiseuille 
flow, which are normalized by the their maximum. 

showing that viscous force term is also proportional to the streamwise pressure gradient. 
The axial velocity is expressed as by integration on above equation, 

Uz =Uq(^- , (17) 

where uq = is the center line velocity, z is in axial direction and r is in radial 

direction of the cylindrical coordinates, and R is the radius of the pipe. The energy gradient 
can be expressed as below for any position in the flow field (noticing Ur = 0), 

Similar to plane Poiseuille flow, the viscous force term is linear, while the kinematic 
energy gradient increases quadratically with the pressure gradient. 

For pipe Poiseuille flows, the ratio of the energy gradient to the viscous force term, K, 
is {dp/dr = 0), 

= d-r[2'^ )^^[-d^^r^)=-^'''['-R-^)/[-^l^ 

pUR r/r^\lr/r^\ , ^ 

' 1 - 7^2 = o^e- 1 - ^ . (19) 



p R\ R^ J 2 R\ 

Here, Re = pUD/ fi, and U = ^uq has been used for pipe Poiseuille flow, uq is the maximum 
velocity at centerline and U is the averaged velocity. 
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Poiseuille pipe 






He = pU U/ fi 




Keynoias (ioooj 




zoUO 




Fetai dz Head (1959) 




zOUO 




Uarbysnire &; Mumn^iyyoj 




1760 ~ z3U0 




Most literature cited 




onnn 
zUUU 


Poiseuille plane 




Re = pUL/fi 


Re = puoh/p 




Davies k White (1928) 


1440 


1080 




Patel & Head (1969) 


1350 


1012 




Carlson et al (1982) 


1340 


1005 




Alavyoon et al (1986) 


1466 


1100 




Most literature cited 




1000 



Table 2: Collected experimental data of the critical Reynolds number for plane Poiseuille 
flow and pipe Poiseuille flow. U is the averaged velocity. uo=2U, and D is the diameter of 
the pipe for pipe Poiseuille flow. uo=1.5U, L =2h, and L is the hight of the channel for plane 
Poiseuille flow. 

The distribution of K along the transversal direction for pipe Poiseuille flow is the same 
as that for plane Poiseuille flow if it is normalized by its maximum and y/h is replaced by 
r/R, and the maximum of K also occurs at r/R = 0.5774, as shown in Fig.7. 

5 Comparison with Experiments 

Experiments for Poiseuille flows indicated that when the Reynolds number is below a crit- 
ical value, the flow is laminar regardless of the disturbances. For circular Poiseuille flow 
(Hagen-Poiseuille) , Reynolds (1883) [2] carried out the first systematic experiment on the 
flow transition and found that the critical Reynolds number for transition to turbulence is 
about Rec = pUD/ji = 2300, where the U is the averaged velocity and D is the diame- 
ter of the pipe. Now, the most accepted critical value is Rcc = pUD/ = 2000 which is 
demonstrated by numerous experiments [28]. All the collected data could be put in a range 
of 1760 to 2300 [9]. There are also experimental data for the transition for plane Poiseuille 
flows in the literature. Davies and White [29] showed that the critical Reynolds number for 
transition to turbulence is Rec = pUL/ p = 1440 for plane Poiseuille flow, where the U is the 
averaged velocity and L = 2/i is the width of the channel. Patel and Head [30] obtained a 
critical value for turbulence transition, Rec = 2000 for pipe Poiseuille flow, and Rec = 1350 
for channel flow through detailed measurements. Carlson et al. [31] found the transition at 
about Rec = 1340 for plane Poiseuille flow using flow visualization technique. Alavyoon et 
al.'s [32] experiments show that the transition to turbulence for plane Poiseuille flow occurs 
around Rec = puoh/p = 1100. The most accepted value of minimum Rec for plane Poiseuille 
flow is about Rcc = puoh/p w 1000 [6] [7]. All the collected experimental data are listed in 
Table 2. Although these experiments are done at various different environmental conditions, 
they are all near a common accepted value of critical Reynolds number. In the following, 
we show that there is a critical value of i^max at which the flow becomes turbulent. In order 
to more exactly compare plane Poiseuille flow to pipe Poiseuille flow at same experimental 
conditions, we prefer here to use Patel and Head's data [30] to evaluate the parameters at 
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Flow type 


Re 


expression 


Eigenvalue analysis 


Energy 


Experiments 


Kmax at Exp 








Rec 


method Rec 


RCc 


Rcc value 


Poiseuille pipe 


Re 


= pUD/fi 


stable for all Re[7] 


81.5 


2000 [30] 


385 


Poiseuille plane 


Re 


= pUL/fi 


7696[10] 


68.7 


1350[30] 


389 




Re 


= puoh/fj, 


5772[10] 


49.6 


1012[30] 


389 


Plane Couette 


Re 


= puhh/n 


stable for all Re[7] 


20.7 


370[37, 38] 


370 



Table 3: Comparison of the critical Reynolds number and the ratio of the energy gradient to 
the viscous force term, K^ax for plane Poiseuille flow and pipe Poiseuille flow. The critical 
Reynolds number by energy method is taken from Schmid and Henningson (2001). 



the critical conditions. Patel and Head's data are also the best to fit all of the data and are 
cited by most literature. 

Now, we calculate the critical value of i^max at the transition condition for both plane 
Poiseuille flow and pipe Poiseuille flow using Eqs.(15) and (19), respectively. For plane 
Poiseuille flow, one obtains -Rrmax= 389 at the critical Reynolds number Rec= 1350. For 
pipe Poiseuille flow, one obtains Kmax= 385 at the critical Reynolds number Rec= 2000. 
These results are shown in Table 3. In this table, the critical Reynolds number obtained 
from energy method is also listed. From the comparison of critical values of -R'max for 
plane Poiseuille flow and pipe Poiseuille flow, we find that although the critical Reynolds 
number is different for the two flows, the turbulence transition takes place at the same -fCmax 
value, about 385 ^ 389. This demonstrated that i^max is really a dominating parameter 
for the transition, and -RTmax is a better expression than the Re number for the transition 
condition. We can further conclude that energy gradient theory is better than the linear 
stability theory for the prediction of critical Reynolds number of subcritical transition. In 
this way, the proposed idea is verified for wall bounded parallel shear flows. Therefore, it 
may be presumed that the transition of turbulence in other complicated shear flows would 
also depend on the i^max in the flow field. 

Nishioka et al (1975)'s famous experiments [11] for plane Poiseuille flow showed details 
of the outline and process of the flow breakdown. The measured instantaneous velocity 
distributions suggest that the break down of the flow is a local phenomenon, at least in 
its initial stage. As in Fig. 8, the base flow is laminar and the instantaneous distribution 
of the velocity breaks at the position y/h = 0.50 (T = 4 to 6) to 0.62 (T = 8 to 9) by 
showing an oscillation of velocity in y/h = 0.50 ~ 0.62. They show an inflectional velocity 
in this range of y/h. This result means that the flow breakdown first occurs in the range 
of y/h = 0.50 ~ 0.62. This coincides to the prediction of our theory, i.e., the position of 
-f^max is the most dangerous point which occurs at y/h = 0.5774. These results are enough 
to confirm the theory of "energy gradient" valid at least for Poiseuille flows (pressure driving 
flow). 

For pipe flow, in a recent study [33], Wedin and Kerswell showed that there is the presence 
of the "shoulder" in the velocity profile at about r/R = 0.6 from their solution of travelling 
waves. They suggested that this corresponds to where the fast streaks of traveling waves 
reach from the wall. It can be construed that this kind of velocity profile as obtained by 
simulation is similar to that of Nishioka ct al's experiments for channel flows. The location 
of the "shoulder" is about same as that for .fCmax- According to the present theory, this 
"shoulder" may then be intricately related to the distribution of energy gradient. The 
solution of traveling waves has been confirmed by experiments more recently [34] . 
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Figure 8: Instantaneous velocity distributions in a plane Poiseuille flow (Nishioka et al. 1975) 
[11]. Time T corresponding to each distribution is noted on the trace of the u fluctuation at 
y/h=0.6, sketched in the figure. Solid circle is the mean velocity. Uc is the velocity on the 
channel center-plane. y/h=0 is at the center-plane and y/h=l is at the wall. Instantaneous 
velocity U+u at a point is composed of the mean velocity U and the fluctuation velocity u. 
(Courtesy of Nishioka; Use permission by Cambridge University Press). 



15 



250 



500 750 1000 R^^^ PoiseuiiiePfane 



(a) 



pUD 



Poiseuille Pipe Re — 
J \ I M 



D » soo 7SQ 1000 mm ism nwa mm 




Q: 250 500 750 1000 1250 1500 1750 2OO0 

Potseujile Pipe 



Figure 9: Scaling of plane Poiseuille flow with pipe Poiseuille flow, (a) Keeping the Re is 
constant, (b) Keeping the Xmax is constant. The correlation should be carried out at the 
same Kmax vaue (b), rather at the same Re (a). 

The energy gradient theory can also be used to explain the reason why it is not appro- 
priate to scale the outer flow or overlap profiles of channel flow and those of pipe flow at the 

same Reynolds number turbulent flows[35]. This is easily understood by the fact that pipe 
Poiseuille flow has a same velocity and energy gradient distributions in the radial direction 
as the plane Poiseuille flow has in the y direction, but the former has a smaller hydraulic 
diameter than the latter and has more viscous friction. Therefore, the scaling should be 
carried out at the same i^'max value, but not at the same Reynolds number (Fig. 8). At a 
same Reynolds number, for example, say, 1000, the plane Poiseuille flow reaches the critical 
Re number for transition, while pipe Poiseuille flow is far from the critical Re number, as 
shown in Fig.Qa. The flow state at these two flows are definitely different at this Re number. 
This principle also applies to turbulence flow range. If we compare the two type of flows at 
same i^max value, they should have the same flow behaviour (Fig. 9b). 

In a separating paper [36] (owing to the space limit here), we apply the energy gradient 
theory to the shear driving flows, and show that this theory is also correct for plane Couette 
flow. We obtain -R'max = 370 at the critical transition condition determined by experiments 
below which no turbulence occurs (see Table 3). This value is near the value for Poiseuille 
flows, 385 ^ 389. The minute difference in the number is not important because there is 
some difference in the determination of the critical condition. For example, the judgement 
of transition is from the chart of drag coefficient in Patel and Head [30], while visualization 
method is used in [37, 38]. These results demonstrate that the critical value of i^max at 
subcritical transition for wall bounded parallel flows including both pressure driven and 
shear driven flows is about 370 ^ 389. 

More recently, the energy gradient theory is applied to the Taylor-Coucttc flow between 
concentric rotating cylinders [39]. The detailed derivation for the calculation of the energy 
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Figure 10: Comparison of the theory with the experimental data for the instabihty condition 
of Taylor-Couette flow between concentric rotating cyhnders (i2i=3.80cm, i22=4.035 cm). 

radius of the inner cyhnder; R2: radius of the outer cyhnder. wi and u;2 are the angular 
velocities of the inner and outer cyhnders, respectively. The critical value of the energy 
gradient parameter Kc=139 is determined by the experimental data at u;2 = and oji ^ 
(the outer cylinder is fixed, the inner cylinder is rotating). With i^c=139, the critical value 
of oJi/v versus 0J2/V is calculated using the energy gradient thoery for UJ2/1' = -2200 — 900 
[39]. 
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gradient parameter is provided in the study. The theoretical results for the critical condition 
of primary instability obtain very good agreement with Taylor's experiments (1923) [40] and 
others, see Fig. 10. Taylor (1923) used mathematical theory and linear stability analysis and 
showed that linear theory obtained agreement with the experiments. However, as is well 
know and discussed before, linear theory is failed for wall bounded parallel flows. As shown 
in this paper, the present theory is valid for all of these said flows. Therefore, it is concluded 
that the energy gradient theory is a universal theory for the flow instability and turbulent 
transition which is valid for both pressure and shear driven flows in both parallel flow and 
rotating flow configurations. 

The Rayleigh-Benard convcctive instability and the stratified fiow instability could also 
be considered as being produced by the energy gradient transverse to the flow (thermal or 
gravitational energy) . The energy gradient theory can be not only used to predict the gener- 
ation of turbulence, but also it may be applied to the area of catastrophic event predictions, 
such as weather forecast, earthquakes, landslides, mountain coast, snow avalanche, motion 
of mantle, and movement of sand piles in desert, etc. The breakdown of these mechanical 
systems can be universally described in detail using this theory. In a material system, when 
the maximum of energy gradient in some direction is greater than a critical value for given 
material properties, the system will be unstable. If there is a disturbance input to this sys- 
tem, the energy gradient may amplify the disturbance and lead to the system breakdown. 
This problem will be further addressed in future study. 

6 Conclusions 

The mechanism for the flow instability and turbulent transition in parallel shear flows is 
studied in this paper. The energy gradient theory is proposed for the flow instability. The 
theory is applied to plane channel fiow and Hagen-Poiseuille fiow. The main conclusions of 
this study are as follows: 

1. A mechanism of fiow instability and turbulence transition is presented for parallel 
shear fiows. The theory of energy gradient is proposed to explain the mechanism of flow 
instability and turbulence transition. It is stated that the energy gradient in transverse 
direction tries to amplify the small disturbance, while viscous friction in streamwise 
direction could resist or absorb this small disturbance. Initiation of instability depends 
upon the two roles for given initial disturbance. Viscosity mainly plays a stable role to 
the initiation of flow instability by affecting the base fiow. 

2. A universal criterion for the fiow instability initiation has been formulated for wall 
shear fiows. A new dimensionless parameter characterizing the flow instability, K, 
which is deflned as the ratio of the energy gradients in transverse direction and that in 
streamwise direction, is proposed for wall bounded shear flows. The most dangerous 
position in the flow field can be represented by the maximum of K. The initiation 
of flow breakdown should take place at this position first. This idea is confirmed by 
Nishioka et al.'s experiments. 

3. The concept of energy angle is proposed for flow instability. This concept helps to 
understand the mechanism of viscous instabilities. Using the concept of energy gradient 
and energy angle, it is theoretically demonstrated for the first time that viscous flow 
with a velocity inflection is unstable. 
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4. It is demonstrated that there is a critical value of the parameter -fCmax at which the 
flow transits to turbulence for both plane Poiseuille flow and pipe Poiseuille flow, 
below which no turbulence exists. This value is about i^max= 385 ^ 389. Although 
the critical Reynolds number is different for the two flows, the turbulence transition 
takes place at the same ifmax value. 

5. The energy gradient theory is a universal theory for the flow instability and turbulent 
transition which is valid for both pressure and shear driven flows in both parallel flow 
and rotating flow configurations. 
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